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ABSTRACT / 

Theoretically determined frequencies and absorption intensities are reported for 
the vibrational spectrum of the covalent HOOOH and hydrogen bonded HO—HOO 
intermediates that may form in the reaction of the hydroxyl and hydroperoxyl radicals. 
Basis sets of DZP quality, augmented by diffuse and second sets of polarization functions 
have been used with CASSCF wave functions. The calculated harmonic vibrational 
frequencies of HOOOH have been corrected with empirical factors and presented in the 
form of a "stick" spectrum. The oxygen backbone vibrations, predicted to occur at 519, 
760, and 870 cm' 1 , are well separated from most interferences, and may be the most 
useful for the species’ identification. In the case of the hydrogen bonded isomer, 
emphasis has been placed upon prediction of the shifts in the intramolecular vibrational 
frequencies that take place upon formation of the complex. In particular, the HO stretch 
and HOO bend of HO z are predicted to have shifts of -59 and 53 cm' 1 , respectively, 
which should facilitate their identification. It is also noted that the antisymmetric 
stretching frequency of the oxygen backbone in HOOOH exhibits a strong sensitivity to 
the degree of electron correlation, such as has been previously observed for the same 


mode in ozone. 
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I. INTRODUCTION 

The reaction of the hydroxyl and hydroperoxyl radicals to produce water and 

oxygen 

HO- + H0 2 - = H 2 0 + 0 2 (I) 

has received a great deal of experimental and theoretical attention in recent years. In its 
own right, this reaction plays an important role in determining hydroxyl radical 
concentrations in models of atmospheric and combustion chemistry. It is also of interest 
as a prototype for reactions involving more complex radical species. The bimolecular 
rate constant for reaction (1) was found to be quite large. Early measurements of it 
exhibited an apparent pressure dependence, with k t increasing from about 
0.1-0.5 x 10' 10 cm 3 s* 1 molecule' 1: 1 ’ 2,3,4 to about 

1-2 x 10' 10 cm 3 s‘ 1 molecule' 1: 5,6,7,8,9,10 as the total pressure was increased from a few 
torr to one atmosphere. The initial investigation of the temperature dependence 11 of k x 
yielded at estimate of E a /R = -416 ± 86 K. The large value of kj coupled with its 
apparent pressure dependence and negative activation energy generated a great deal of 
interest 12,13,14,15 in the detailed mechanism of reaction (1), with considerable debate 
over whether it followed a direct abstraction mechanism (possibly via a hydrogen bonded 
intermediate), 

HO- + H0 2 - ~ H0-H0 2 - H0 -H--0 2 - H 2 0 + 0 2 (la) 

proceeded via a stabilizable covalent intermediate, 
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h 2 o + o 2 W 

HO- + H0 2 - ** [HOOOH]* 

^ M 

H 2 0 3 (ic) 

or followed both pathways. 

Recent studies 16,17 have indicated that the low and high pressure measurements 
of kj may agree to within the precision of the data, and the apparent pressure 
dependence is currently explained in terms of experimental uncertainties and systematic 
errors in the earlier data. The temperature dependence of reaction (1) has also been 
reinvestigated 16 , indicating a smaller value of E a /R=-250 ± 50 K. The rate constant 
value at 298 K currently recommended by the NASA panel 18 for use in atmospheric 
modeling (1.1 x lO'^cmV'molecule' 1 ) has no pressure dependence and a temperature 
dependence given as E a /R = -250 ± 200 K. 

Two isotopic substitution studies have been carried out to elucidate the 
mechanism of reaction (1). Kurylo et al ? studied the H ls O + H0 2 ’ reaction at high 
pressure (1 atm of SF 6 ) and concluded that formation of the products did not proceed 
via the HOOOH intermediate of reaction (lb), because an insufficient amount of the 
16 0 18 0 product was found. A low pressure H ls O study by Dransfeld and Wagner 13 
found the rate for the reaction 

H k O- + H0 2 - H 2 ls O + 0 2 
to be roughly equal to that of the exchange reaction 

H ls O- + H0 2 - H l6 0- + H m 0 l6 0-. 

These authors attributed the formation of water and molecular oxygen to mechanism 
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(la) and indicated that the formation of the HOOOH intermediate resulted in 
dissociation to reactants rather than stabilization under their conditions. The results of 
Kurylo et al? were reconciled with their own by suggesting that the stabilization path (lc) 
completely dominated the isotope exchange dissociation under the conditions of Ref.7 (1 
atm of SF 6 ), thus accounting for the fact that very little 16 0 18 0 had been observed. 

These mechanistic studies, combined with theoretical results and the knowledge 
that the pressure dependence is probably not strong, has encouraged analysis of the 
mechanism of reaction (1) in terms of the direct abstraction pathway. 15,17,19,20 The 
existence of a potential energy well corresponding to a weakly bound hydrogen bonded 
intermediate, as in reaction (la), can be used to account for the negative temperature 
dependence of reaction (1). 

Early theoretical work on this system focused on the ground state potential energy 
surface of the HOOOH intermediate. Work from this laboratory 21 and others 15,19,20 has 
provided stability estimates of both the covalently bonded and hydrogen bonded 
intermediates. The earliest work 22 that included electron correlation paid particular 
attention to the potential surface for torsional motion about the 0-0 single bonds in 
HOOOH. Recent computational studies 15,19,20 have characterized not only the energy 
minima of the intermediates, but also transition states on some of the low lying surfaces. 
Gonzalez et al. 19,20 have employed ab initio potential surfaces ( ! A, 3 A’, and 3 A") in 
combination with rate theory calculations to study both the temperature and pressure 
dependence of reaction (1). Their results generally are in agreement with the 
experimental interpretations: product formation at normal conditions is due almost 
entirely to reaction (la) on the triplet surfaces; the temperature dependence of reaction 
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(la) can be modeled in reasonable agreement with experiment; the pressure dependence 
due to reaction (lc) will not be important with helium as bath gas below about 5 
atmospheres; and, due to a significant activation energy, rearrangement of HOOOH on 
the ! A surface to form the overall products, reaction (lb), will not be important except at 
high temperatures. 

It is expected that reaction (1) will continue to receive both theoretical and 
experimental attention and that the nature of the intermediates will remain central to 
these investigations. While the present evidence indicates 19 that in helium the HOOOH 
intermediate is not stabilized except at pressures exceeding five atmospheres, the 
suggestion 13 stabilization is dominant over unreactive dissociation in one atmosphere of 
SF 6 raises questions about the importance of channel (lc) in atmospheric and 
combustion environments. In comparing the ratio of k u to the rate constant for 
formation of [HOOOH]’, the authors of Ref. 20 find a 6-fold disagreement between 
value deduced from the isotope studies and that arising from the rate theory calculations. 
In further studies at high pressure it would be very interesting to probe directly for the 
presence of stabilized intermediates; the high pressure flow system recently developed by 
the Harvard group 23 is especially intriguing in this regard. Since Fourier transform 
infrared spectroscopy (FUR) could be an effective detection method for any such 
intermediate, it would be useful to have in hand theoretical predictions of the vibrational 
absorption spectrum (frequencies and intensities) for these intermediate species. 

From a theoretical point of view, the vibrational frequencies of HOOOH are 
interesting because of the molecule’s similarity to ozone. The magnitude and order of 
the calculated ozone vibrational frequencies is extremely sensitive to the level of theory 
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employed. 24 In that case the well known requirement to describe adequately the 
biradical character of the ground state 25 becomes entangled with other aspects of the 
electron correlation problem, and a clear understanding of the difficulty in describing the 
vibrational frequencies is not easily obtained. The l A ground state of HOOOH has no 
biradical complications and is reasonably well described by a closed shell wave function, 
making examination of its -OOO- backbone modes with those of ozone quite interesting. 

The present study reports the results of an ab initio quantum chemical 
determination of the vibrational frequencies and absorption intensities of both the 
covalently bonded HOOOH intermediate and the triplet ground state hydrogen-bonded 
HO— H0 2 species. Vibrational frequencies calculated using many-body perturbation 
theory methods have appeared in the literature for HOOOH, but the intensities have not 
been previously reported. Neither the vibrational frequencies nor the absorption 
intensities have been reported for the 3 A’ ground state of the hydrogen-bonded complex, 
although the frequencies of the closely related 3 A" excited state minimum and transition 
state have been published. 20 
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II. COMPUTATIONAL METHOD 

A. Basis Sets. 

This study employed a polarized double-zeta (DZP) quality Cartesian Gaussian 
basis set consisting of Dunning’s [4s2p 1 2s] contraction 26 of Huzinaga’s (9s5p 1 4s) set. 27 
A scale factor of C = 1.2 was chosen for the hydrogen s functions, and polarization was 
provided by a set of d functions on oxygen (exponent, a =0.85) and a set of p functions 
on hydrogen (exponent, a = 1.0). In some phases of the work the DZP basis was 
augmented by replacing the single set of d functions with a double set (exponents, 
a j = 0.425, a 2 = 1.70) and is denoted DZP(2d). When an sp set of diffuse functions was 
present on the oxygen atoms (exponent, a =0.059), the basis set notation is given as 
DZP( + ) or DZP(2d, + ). Because of program limitations, some of the many-body 
perturbation theory calculations were carried out with a split valence 6-3 1G’ basis 28 
that is of similar quality to the DZP basis. 

B. Electronic Structure Calculations. 

In general the SCF calculations reported here are of either closed- or open-shell 
spin restricted type, except that UHF calculations were utilized to define the orbital basis 
for the open-shell perturbation theory results. The MCSCF calculations used throughout 
this study are of the complete active space 29 (CASSCF) variety and utilize a density 
matrix driven algorithm. 30 The notation "MC[n,m]" will be used to designate an 
expansion that includes all possible configurations for n electrons in m molecular 
orbitals. Geometry optimization was performed using the algorithm of Baker. 31 The 
many-body perturbation theory (MBPT) calculations were carried out at the second and 
fourth order (MP2 and MP4[SDTQ]) using the methodology of Pople and coworkers. 32 
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C. Vibrational Spectra. 

The vibrational spectra were calculated in the "double harmonic" 
approximation 33 in which the potential surface is assumed to be quadratic and the 
electric dipole moment is assumed to vary linearly near the potential surface minimum. 

In most cases the second derivatives of the energy were calculated by numerical central 
differences of analytical energy gradients that had been determined at small 
displacements about equilibrium. In the case of SCF wave functions, the energy second 
derivatives were directly calculated by analytical methods, and in the MP4 case they were 
determined as second differences of energies at displaced geometries. 

D. Computer Programs. 

The MCSCF and SCF calculations were carried out using the North Dakota State 
version of the GAMESS 34 electronic structure package, and the MBPT calculations 
utilized the Gaussian 86/88/90 programs. 35 Computers used in this study included the 
Convex Cl, the Sun Sparc 2, and the Cray YMP. 
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III. RESULTS AND DISCUSSION 
A. The Covalent HOOOH Species. 

1. Equilibrium geometry 

The results of this study for the minimum energy conformation of the *A ground 
state of the HOOOH molecule are reported in Table I, along with those previous results 
that employed correlated wave functions. The molecule displays C 2 symmetry at the 
energy minimum, and from the table it is noted that the bond lengths and angles are 

r 



HOOOH Covalent Intermediate 

quite insensitive both to the level of the CAS treatment and to augmentation of the basis 
set. The active space of the MC[8,8] wave functions includes the eight sigma bonding 
electrons in the space of the four bonding and four antibonding molecular orbitals. The 
MC[14,11] wave functions include one nonbonding lone pair on each oxygen in the active 
space, giving an MCSCF expansion of 32,670 configuration spin functions at unsymmetric 
(Cj) conformations. 

In general with basis sets of this quality, it is expected that experimental bond 
length would be bracketed by the SCF and CASSCF level results, and the experimental 
value of R oh should lie in the range 0.965-0.970 A. However, the OO bond distance 
merits special attention because the change from the RHF to CASSCF results is 
unusually large ( + 0.09 A). To assist in evaluation of these results, we carried out 
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calculations using the same basis set and comparable electron correlation treatments on 
hydrogen peroxide, which has a similar 0-0 single bond. In that case, the optimum 0-0 
bond length was found to increase from 1.391 to 1.487 A as the treatment progressed 
from RHF to MC[ 10,8], which is comparable to the progression from RHF to MC[14,11] 
for H 2 0 3 . This R 00 value of 1.487 A is about 0.03 A longer than the most accurate H 2 0 2 
theoretical 36 and experimental 37 values, 1.462 and 1.464 A, respectively. By analogy, it 
is expected that a more accurate treatment of H 2 0 3 including dynamic electron 
correlation would reduce the value of R 00 to about 1.43 A, which is predicted to be 
quite close to the experimental value. 

Based upon similar calculations for H 2 0 2 and H z O, the HOO bond angle from 
our MC[14,11] calculations are expected to be reasonably accurate, with the 
experimental result predicted to be several tenths of degree larger. The OOO angle is 
expected to be sensitive to the correlation treatment of the oxygen lone pairs, and the 
results in Table I show that it varies considerably as the MCSCF active space is enlarged. 
However, the overall variation in 0 OOO is only 2.5°, so we predict that the MC[14,11] 
values are probably accurate to within ± 1°. The dihedral angle is harder to predict 
because of the flatness of the torsional potential and may be in error by several degrees. 
2 Vibrational frequencies 

The harmonic vibrational frequencies calculated for the covalent H 2 0 3 species are 
presented in Table II, along with the available experimental observations 38 . The 
subscripts used to index the normal modes do not follow the sequence in Ref. 38. 

Rather, they are assigned by placing the normal modes in order of increasing frequency 
as determined at the RHF/DZP level. Examination of the normal modes provided the 
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following descriptions: v j (symmetric torsional rotation); v 2 (antisymmetric torsional 
rotation); v 3 (OOO bend); v 4 (symmetric 00 stretch); v 5 (antisymmetric OO stretch); 
v 6 (symmetric HOO bend); v 7 (antisymmetric HOO bend); v 8 (antisymmetric OH 
stretch); and v 9 (symmetric OH stretch). Within each column of Table II are the 
frequencies for a particular normal mode as calculated with a variety of basis sets and 
wave functions. Thus, if the frequencies of two normal modes switch order as the 
computational level is changed, some rows of the table will not be in strictly increasing 
order. 

Use of the DZP(2d) basis set was not extended to the rather costly MC[14,11] 
calculations, because comparison of DZP and DZP(2d) frequencies at the RHF and 
MC[8,8] levels had already shown (Table II) that the extra set of d functions had an 
effects of only about 1%. There was no reason to anticipate a larger basis set effect 
using the MC[14,11] wave function. The MP2 frequencies were taken from the 
literature, but the MP4 results were added by us for comparison. 

The torsional frequencies v : and v 2 decrease by about 15% in going from RHF to 
MC[8,8] wave functions and then increase by as much as the same amount when the 
oxygen lone pairs are included in the active space (MC[14,11]). These torsions are 
expected to be highly very anharmonic, so that their predicted frequencies are not 
expected to be very accurate. However, the large changes between MC[8,8] and 
MC[14,11] results indicate that an accurate description of these torsional motions about 
the 0-0 bonds should include a correlated description of the oxygen lone pairs. 

In the cases of "normal" bending and stretching vibrations, harmonic frequencies 
calculated at the DZP/RHF level frequently overestimate the observed fundamentals by 
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10-15%, with the overestimation being reduced to just a few per cent with correlated 
treatments. This residual overestimation is due to neglect of the anharmonic force 
constants and basis set limitations, with error cancellations often yielding errors 
somewhat smaller than might be expected within the harmonic approximation. As shown 
in Table II, modes 6-9, which involve hydrogen atom motions, follow the usual pattern 
with the expected decreases in frequency in progressing from RHF to MC[8,8] 
treatments, and only small changes as the active space is increased to MC[14,11], 

The oxygen backbone vibrations, normal modes 3-5, show behavior that departs 
from the "usual" pattern, with relatively large frequency changes being noted as the 
computational level is varied. For example, Table II shows that with the DZP basis set, 
the decreases in frequency from RHF to MC[8,8] are quite large: 19% for v 3 , 29% for 
v 4 , and 31% for v 5 . In further expanding the treatment of electron correlation the 
MC[ 14,11] wave function, the OOO bending motion (v 3 ) is unaffected, the symmetric 
OO stretching frequency (v 4 ) increases by about 3%, and the antisymmetric OO 
stretching frequency (v 5 ) further decreases by about 11%. In going from MC[8,8] to 
MC[14,11] wave functions, the oxygen stretching frequencies not only exhibit rather large 
changes, but also switch order. While both MPBT treatments place the frequencies in the 
same order as the larger CAS calculations (v 4 >v 5 ), the changes in going from MP2 to 
MP4 are unusually large (3-6%) for modes 3-5, again indicating a high sensitivity to the 
treatment of the correlation energy. 

Examination of the MC[14,11] natural orbitals shows that the three nonbonding 
electron pairs being included in the active space are best described as a 2s-type on the 
central oxygen atom and a pair of 2p-type on the noncentral oxygens. Attempts to 
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replace the central oxygen 2s-lone pair with its 2p-lone pair either failed or converged to 
a higher-energy solution. 

3. Comparison with hydrogen peroxide 

One tool for critical evaluation of the data presented in Table II involves 
examination of similar computational results for related systems that have well known 
experimental spectra. In Table III are displayed results of calculations we have carried 
out for hydrogen peroxide using the same basis sets and comparable MCSCF active 
spaces. For H 2 0 2 the MC[6,6] and MC[10,8] active spaces are analogous, respectively, to 
the MC[8,8] and MC[14,11] spaces of H 2 0 3 . The smaller active space includes only the 
six bonding electrons in the three sets of nominally bonding and antibonding orbitals, 
whereas the larger one includes one nonbonding lone pair on each oxygen atom. Table 
III also shows results for H 2 0 2 labelled ”MC[14,10]", which include in the active space 
both lone pairs on each oxygen atom, a calculation which we were not able to carry out 
on the larger H 2 O a system. The theoretical H 2 0 2 harmonic frequencies and comparable 
MP2 values 39 are reported along with their ratios to the observed fundamentals. 40 The 
MP2 results were chosen for inclusion because they are comparable in both basis set 
quality and neglect of anharmonicity to the present work; there are much better 
theoretical values for H 2 0 2 available from the same authors and others. 36,39 

From the ratios in Table III it is seen that the frequencies of the H 2 0 2 normal 
modes that primarily involve hydrogen motion generally follow the usual pattern, with 
10-12% decreases in going from RHF to the MC[6,6] level and only a very small change 
as the active space is further expanded. The frequency changes reported in Table II for 
v 6 -v 9 , the HOO bends and HO stretches in H 2 0 3 , follow very much this same pattern. It 
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is expected that their residual errors will be similar to those for H 2 0 2 . The 0-0 stretch 
in hydrogen peroxide shown in Table III is overestimated by 33% at the RHF level and 
underestimated by 5-6% at each of the MCSCF levels. The large change in frequency 
from RHF to MC[6,6] mirrors closely that of H 2 0 3 oxygen-oxygen stretches, v 4 and v s in 
Table II. However, the small subsequent change in the H 2 0 2 0-0 stretch (Table III) in 
going from MC[6,6] to MC[10,8] is in sharp contrast to the situation for v 4 and v 5 of 
H 2 0 3 (Table II) in proceeding from MC[8,8] to MC[14,11]. The large changes and 
reordering of the OO stretching frequencies of H 2 0 3 upon correlation of the oxygen lone 
pairs do not seem to have precedents in the H 2 0 2 results. 

4. Comparison with ozone 

To further assist in evaluation of the calculated H 2 0 3 oxygen backbone vibrational 
frequencies, it is instructive to consider the vibrational frequencies of ozone, which are 
well known to be difficult to predict. Lee and Scuseria 24 have recently reviewed the 
literature on this topic. The ground state of ozone has significant biradical character 25 
and requires a wave function of at least two configurations for a qualitatively correct 
description. The experimentally derived 41 harmonic vibrational frequencies of ozone 
are 1135 cm' 1 (symm. str.), 716 cm' 1 (bend), and 1089 cm" 1 (antisymm. str.). At the 
DZP/SCF level, calculated harmonic vibrational frequencies are found 42 to be much 
too large: 36, 18, and 33%, respectively. Improvement of the wave function to the two- 
configuration level 42 (TCSCF), which provides a qualitatively correct description of the 
molecule’s biradical character, brings about considerable improvement in the bend and 
symmetric stretch frequencies (now 9% and 5% too large), but worse agreement for the 
antisymmetric stretch (38% too large), placing the two OO stretches in the wrong order. 
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Systematic enlargement of the basis set to the TZ2P level does not result in any 
significant improvement with the TCSCF wave function. Inclusion of external electron 
correlation via either configuration interaction (CISD) based upon both the one- or two- 
configuration reference states 43 or perturbation theory 44 (MP2 and MP4) based upon a 
single configuration reference did not result in significant improvement and sometimes 
yielded results clearly inferior to the uncorrelated ones. 

The only approaches to the ozone problem that have yielded significantly 
improved results are those that include a fairly powerful treatment of the nondynamical 
electron correlation. A DZP/CASSCF wave function 45 with the Is and 2s electrons 
inactive gave anharmonic stretching frequencies that were both too small, but placed in 
the correct order. The relative errors in comparing the calculated anharmonic 
frequencies with the observed fundamentals are: -5% (symm. str.), -2.5% (bend), and - 
11% (antisymm. str.). In another DZP/CASSCF study 46 with only the Is electrons 
inactive, a trajectory method was used to determine the anharmonic vibrational 
frequencies for the two stretching motions with the OOO angle fixed at its equilibrium 
value. When compared to the observed fundamentals, these anharmonic frequencies 
were found to be in the correct order and in error by only -4% (symm. str.) and -2% 
(antisymm. str.), a substantial improvement for the antisymmetric stretch. Several 
coupled-cluster based approaches have been applied to the ozone problem. By 
incorporating triple excitations via the CCSD(T) method and using a large ANO basis set 
including /-functions, Lee and Scuseria 24 were able to predict the harmonic frequencies 
of ozone in the correct order and with errors of less than 3%. 

Similarities are apparent between the data in Table II for modes v^, the oxygen 
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backbone vibrations of H 2 0 3 , and the 0 3 vibrations. Regardless of basis set quality, the 
frequencies of modes 3-5 of H 2 0 3 decrease much more in going from RHF to correlated 
levels, than the usually expected ~ 10%. As the wave function is improved from the 
MC[8,8] to MC[14,11] level, which enlarges the active space to include one lone pair on 
each oxygen atom, the antisymmetric stretching frequency v 5 drops greatly and switches 
relative order with the symmetric stretch v 4 . This is reminiscent of the antisymmetric 
stretching frequency in ozone, which takes on wildly differing values and also exchanges 
order with the symmetric stretch as the level of the theoretical treatment is varied. Just 
as the DZP/CASSCF treatment of ozone yielded frequencies with modest errors that 
were in the correct order, we expect that the H 2 0 3 MC[14,11] results are dependable. In 
particular, the large changes in v 4 and v 5 between MC[8,8] and MC[14,11] are thought to 
be significant and not artifacts. 

The MP2 and MP4 results which are reported in Table II for comparison, may be 
somewhat too high, but are probably not greatly in error, as they were in the ozone case. 
The single reference configuration upon which the MBPT treatment is based is a much 
better approximation to the H 2 0 3 wave function than it is to the essentially biradical 
ozone ground state. Also presented for comparison are the experimental fundamentals 
which are available from a study 38 of H 2 0 3 trapped in a condensate. Because the 
aqueous condensate would be expected to be hydrogen bonding to the H 2 0 3 , we are 
hesitant to apply these frequencies directly to the predicted gas phase spectrum. 

However, it is noted that these experimental results are not inconsistent with the 
expectation that the DZP( + )/MC[14,ll] level oxygen-oxygen stretching frequencies may 
be in error by only a few percent. 
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5. Absorption intensities 

Calculated H 2 0 3 infrared absorption intensities are tabulated in Table IV. In 
going from RHF to correlated calculations, there are considerable changes, particularly 
for the HO stretches v 8 and v 9 . Part of this change can be accounted for by a 
concomitant decrease in the electric dipole moment. The intensity of mode 5, the 
antisymmetric stretching vibration, undergoes a large decrease in going from RHF to 
MC[8,8] wave functions, but then experiences even larger increases as the calculation is 
further enlarged to MC[14,11]. Just as the frequency of this mode is quite sensitive to 
the treatment of electron correlation, so is its intensity, with the MC[14,11] calculations 
showing it to be the strongest of the oxygen backbone modes. 

The overall picture is of two very intense low frequency torsional vibrations 
(v j and v 2 ), two weak symmetric stretching vibrations(v 4 and v 9 ), and five vibrations of 
moderate intensity. Augmentation of the basis set from DZP to DZP(2d) does not have 
a large effect on the pattern of intensities. However, addition of a diffuse function in the 
DZP( + ) basis set is important at the MC(14,11) level, with the intensity of mode 8 
(antisymm. OH str.) being doubled and that of mode 3 (-OOO- bend) being reduced by 
one-third. The MP2 intensities, calculated for comparison, are in reasonable agreement 
with the MC[14,11] results, 
d Predicted vibrational spectrum 

In order to present a theoretical spectrum with the greatest predictive value, we 
have empirically corrected the calculated frequencies where we can reasonably do so. 

All of our calculated frequencies are obtained in the harmonic approximation, and this 
neglect of the anharmonic aspects of the vibrational potential should, by itself, cause the 
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predicted frequencies to be several percent above the observed fundamentals. However, 
other shortcomings of the calculation, such as the limited treatment of electron 
correlation, may cause the harmonic frequencies to be underestimated, with the 
cancellation of these errors bringing the calculated and observed values into better 
agreement. We will look to the hydrogen peroxide and ozone examples to obtain 
estimates of the residual errors in the calculated frequencies. 

The estimated absorption intensities are probably less accurate than the 
frequencies, since they rely on detailed aspects of the charge distribution as well as on 
the shape of the potential energy surface. As an example, we note the strong variation 
in intensity in going down the v 8 column of Table TV even though the corresponding 
frequency data follow the usual predictable pattern (Table II). General trends in 
absorption intensities are not as well known, and comparison to experimental numbers is 
not that straightforward. Consequently, we will base our predicted absorption spectrum 
on the uncorrected H 2 0 3 intensities as calculated at the DZP( + )/MC[14,ll] level, our 
best overall calculation. These intensities should be useful for prediction of overall line 
strength patterns, but should not be expected to distinguish between lines of similar 
intensity. 

The torsional vibration frequencies of H 2 0 3 , v t and v 2 , are left uncorrected 
beyond the values obtained at the DZP( + )/MC[14,ll] level. These vibrations may be 
expected to have large anharmonic corrections, but will probably not be well 
characterized by the trends usually found for normal stretches and bends, making a 
systematic correction impractical. Cremer 22 has shown that barriers of only a few 
kcal mol' 1 separate local minima along the path of torsional motion, making it likely that 
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the spectrum will contain complications arising both from torsional tunnelling and 
coupling to the other vibrational modes. Detailed treatments such as have been carried 
out for the torsional mode of H^ 36,39 will probably be needed in order to go much 
beyond the harmonic estimates given here. The (uncorrected) values used for v t 
and v 2 are 337 and 391 cm' 1 . 

The HOO bending and HO stretching frequencies, v 6 -v 9 , are expected to be well 
behaved and, therefore, overestimated at the RHF level by about 10-15% and at the 
CASSCF level by only a few percent. The comparison data in Table III shows that for 
hydrogen peroxide the overestimation at the MC[ 10,8] level is 3-4%. Anticipating the 
same type of behavior for HO stretches and HOO bends in H 2 0 3 , we apply a 4% 
reduction to the DZP( + )/MC[14,ll] values in Table II, arriving at the frequencies 1367, 
1376, 3566, and 3563 cm' 1 for v 6 through v 9 . 

While highly accurate harmonic frequencies will generally be higher than the 
observed fundamentals because of the systematic error due to neglect of anharmonicity, 
less accurate harmonic frequencies determined at the DZP/CASSCF level may be a few 
percent above or below them. The H 2 0 2 data in Table III provides an example where 
the calculated 0-0 stretching mode is about 5% below the experimental values. The 
H 2 0 3 oxygen backbone vibrations have already been noted as not fitting the "usual" 
pattern for well behaved stretches and bends, making uncritical application of the usual 
corrections somewhat problematic. It has also been observed by us that the behavior of 
the H 2 0 3 modes with improvement of the wave function is similar to that of the ozone 
vibrations. In particular, the large decrease in frequency from RHF to MCSCF level and 
the continued large decrease in the antisymmetric stretch frequency as the CASSCF 
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active space is enlarged are similar to the effects found for ozone. Accordingly, we 
consider the previous ozone results in order to formulate corrections to modes 3-5. 

In the previous ozone work a CASSCF calculation 45 comparable to the present 
one yielded anharmonic frequencies that underestimate the experimental fundamentals by 
3-5% for the bend and symmetric stretch and by about 11% for the antisymmetric 
stretch. Use of harmonic frequencies in this comparison would have reduced the 
apparent disagreement by several percent due to cancellation of errors. In fact, the 
harmonic frequencies cited 43 as derived from the work of Ref. 45 differ from the 
observed fundamentals by: -0.5% (symm str), -2% (bend), -5% (asymm str). Considering 
both these previously reported errors for ozone and that found by us for the 0-0 stretch 
in hydrogen peroxide (Table III), we make the somewhat ad hoc assumption that the 
H 2 0 3 DZP( + )/MC[14,ll] frequencies of Table II are too low and in error by -2% 

(symm str and bend) and -5% (asymm. str). Based upon these error estimates, the 
empirically corrected values of v to be used in predicting the spectrum are, 
respectively, 519, 870, and 760 cm' 1 . Clearly, this correction procedure is empirical, but 
we believe it provides estimates of the experimental frequencies that are within 2-3% of 
the correct gas phase values. 

Figures 1 and 2 show a predicted H 2 0 3 spectrum in "stick form", with the nine 
vibrational frequencies v r v 9 corrected as described above and the intensities unmodified 
from the DZP( + )/MC[14,ll] results in Table IV. To recognize probable experimental 
interferences, we also include the observed fundamental frequencies taken from the 
experimental literature for several other species (HO* 47 , H0 2 * 48,49,53 , H 2 0 50 , and 
HPj 40 ) that are likely to be in the reaction flow tube during the experiment. We obtain 
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absorption intensities to accompany the experimental frequencies using CASSCF 
calculations that employ basis sets and active spaces comparable to those used for H 2 0 3 . 
Depending upon the experimental conditions, the concentrations of these interfering 
species may greatly exceed that of any H 2 0 3 formed, so that even intrinsically weak lines 
could obscure nearby absorptions due to H 2 0 3 . 

Examination of Figures 1 and 2 shows that the OH stretches and HOO bends of 
H 2 0 3 may well be obscured by other species, depending upon the reaction conditions. 
However, the vibrations of the H 2 0 3 oxygen backbone (500-700 cm' 1 ) should not suffer 
significant interference and should be easier to identify. The torsional vibrations are 
predicted to be very intense. However, their frequencies are probably not well predicted 
by the present harmonic treatment, and it is not clear how close they might be to the low 
frequency absorption of hydrogen peroxide. 

B. The Hydrogen-Bonded HO— HOO Species. 

1. Equilibrium geometry 

In previous studies, 15,20,21 the ground state of the hydrogen bonded complex was 



HO **• HOO Hydrogen-bonded Intermediate 
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found to be planar and display 3 A’ symmetry, with the unpaired electrons on each radical 
situated in out-of-plane a” orbitals. The predictions made by these authors of the 
stability of the complex relative to the separated radicals fell in the range 4.5 - 9.0 kcal 
mol' 1 . The 3 A’ state is treated here even though it does not correlate through planar 
conformations with the ground state products of reaction (1). We are focusing on 
intermediate species that might be stabilized and detected, regardless of whether or not 
they are likely to proceed to products. The MC[10,9] calculations reported for 
HO—HOO in Tables V-VI have active spaces that include ten electrons in three a’ (o- 
type) bonding orbitals, three a’ (a '-type) antibonding orbitals, and three a" out-of-plane 
(it and it’ of H0 2 -; 2p of HO) orbitals. None of the in-plane a’ lone pairs on the 
oxygens are included in the active space. 

The minimum energy conformation of the HO—HOO complex found at the 
various levels of theory in the present study is consistent with the previous work. As 
shown in Table V, the intramolecular coordinates are found to be only slightly perturbed 
from those of the uncomplexed radicals. This is even true for the R 00 coordinate of 
H0 2 , which is highly sensitive to the treatment of electron correlation. The MP2 and 
MCSCF results differ by -0.04 A and bracket the experimental value 51 of 1.331 A. 

Even though the various R 00 values in Table V differ by as much as 0.05 A, the changes 
in R 00 upon complexation fall in the range 0.002-0.008 A for each of the methods 
tabulated. 

As anticipated because the potential well of the complex is quite shallow, its 
intermolecular coordinates display greater sensitivity both to the basis set used and to 
the treatment of electron correlation. These parameters would probably undergo 
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considerable change as further improvements in the basis set and electron correlation 
treatment are made. One interesting difference between the results obtained by the 
different methods is that the Ro... H distance increases in going from the ROHF to the 
CASSCF treatment, whereas it decreases in going from ROHF to MP2. While the 
present study has not considered the 3 A" state of the complex, previous studies 15,20 have 
shown that these intermolecular parameters may be considerably different for that 
electronic state, in which the unpaired electron of the hydroxyl radical is rotated into the 
plane of the complex. 

2 Vibrational frequencies 

In table VI are given the harmonic vibrational frequencies as calculated for the 
hydrogen bonded HO— HOO species. The intermolecular vibrational modes v 1 -v 5 , 
which appear in the spectrum only upon complex formation, are probably not well 
described within the harmonic approximation, and their frequencies may not be 
accurately predicted. While all of the intermolecular mode frequencies increase as 
electron correlation is included in the calculation, modes 2 and 5 show unusually large 
increases. If these calculated results are accurate and modes 2 and 5 appear in the 
range of 700-800 cm' 1 , as indicated by the MC[10,9] calculations, they would be at higher 
frequency than the OOO bend and near the antisymmetric 0-0 stretch in the covalent 
HOOOH species. 

Normal modes 2 and 5 are shown, respectively, in Figures 3 and 4 at both the 
ROHF and MC[ 10,9] levels. These modes are seen to involve mostly out of plane 
motions of the hydrogen atoms. At the ROHF level, mode 2 is mainly a rocking motion 
of the hydroxyl radical (Fig. 3) and mode 5 (Fig. 4) is primarily a rocking motion of the 
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hydroperoxyl radical. At the MCSCF level, these modes show coupled motions of the 
two radicals. Mode 2, while still primarily a rocking motion of the OH radical, includes 
an out-of-phase rock of the H0 2 - species. In the same way mode 5 changes from almost 
pure H0 2 - motion to an in-phase rocking of both radicals. Because use of a correlated 
wave function had such a large effect on the frequencies of modes 2 and 5, the 
DZP/MP2 method was also employed to provide a comparison. While the effect is not 
as large as with the CASSCF approach, modes 2 and 5 still undergo substantial increases 
and mode 4, an in-plane motion, shows an unusually large change. 

It is not clear from examination of the wave functions, charge distributions, and 
force constants why these a” out-of-plane modes are affected so much by inclusion of 
electron correlation. Two of the three leading configurations after the SCF reference 
one, however, involve the n -> n" excitation on H0 2 , either by itself as a single 
excitation or as part of the an - a*n' double excitation. Since the n orbital is localized 
mostly on the central oxygen, and the n ’ orbital mostly on the terminal oxygen, these 
important excitations move a significant amount of negative charge to the terminal 
oxygen of H0 2 -. This increased charge would be expected to enhance the electrostatic 
attraction between the hydroperoxyl terminal oxygen and the hydrogen end of the 
hydroxyl radical, resulting in a larger force constant for out-of-plane motion of that 
hydrogen atom. In fact the out-of-plane force constant for the hydroxyl hydrogen 
increases by about 30-fold, from 0.0006 to 0.018 hartree bohr' 1 in going from ROHF to 
MC[10,9] wave function. However, this explanation is at best incomplete, as there is not 
a similar increase in charge at the hydroxyl oxygen to account for the 4-fold 
enhancement of the out-of-plane force constant for the hydroperoxyl hydrogen atom with 
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inclusion of electron correlation. 

The HO—HOO intramolecular modes have descriptions and frequencies very 
much the same as those of the isolated radicals. In the spectra of weakly interacting 
systems, it is often the shifts of these similar vibrational frequencies upon complex 
formation that are most easily characterized. Experience in this laboratory 52 in 
comparing calculated and experimental frequency shifts for the water dimer, has 
indicated that, although the absolute frequencies are much improved upon inclusion of 
electron correlation, the frequency shifts are more accurately calculated at the SCF level 
than at MCSCF levels comparable to those used in this work. It was found in that study 
that CASSCF results frequently underestimated the shifts by large amounts, sometimes 
assigning them the wrong sign. The MP2 results in that study tended to very badly 
overestimate the frequency shifts, with only the shifts calculated at the SCF level having 
consistently correct signs and reasonable magnitudes. Further, these SCF level frequency 
shifts were not very sensitive to augmentation of the basis set beyond the DZP level. 

In Table VII are tabulated the differences between the harmonic vibrational 
frequencies in the HO—HOO complex and in the isolated radicals, as determined in 
parallel calculations. Following our experience 52 with the water dimer, we expect that 
the ROHF/DZP( + ,2d) shifts are the most accurate predictions in Table VII. The 
frequencies with large shifts, particularly the HO stretch and HOO bend in H0 2 -, are 
likely to have the greatest diagnostic value. 

The large difference between the CASSCF and MP2 values for v 6 (the 0-0 
stretch) reported for the complex in Table VI is similar to the difference found between 
the same methods for isolated H0 2 \ In that case, the experimental fundamental 53 of 
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the 0-0 stretch is 1098 cm' 1 and the CASSCF value (1048 cm' 1 ) is in much better 
agreement. Presumably, the CASSCF value for v 6 in the complex is also closer to the 
experimental one. 

3. Absorption intensities 

In Table VIII are shown the calculated absorption intensities for the hydrogen 
bonded complex of H 2 0 3 . All of the intermolecular vibrations of the complex, modes 1- 
5, have significant intensities, with the two out-of-plane (a") vibrations (modes 2 and 5) 
and the in-plane rotation of the hydroxyl radical (mode 4) being particularly strong. Of 
the intramolecular vibrational modes, the intensities of modes 6-7, the 0-0 stretch and 
HOO bend in H0 2 -, do not change much upon complexation. However, the HO stretch 
of HO z - (mode 8) is predicted at all levels of theory to become 3-4 times more intense 
in the complex than in the isolated radical. The hydroxyl radical HO stretch (mode 9) is 
also predicted by all of the calculations to grow in intensity upon complex formation. 

The experimental intensity of the HO- vibration is weak, and the MCSCF results for the 
isolated radical are probably more accurate than are the ROHF results, suggesting a 15- 
fold increase in intensity upon complex formation. 

The decrease in intensity for mode 2 in going from ROHF to MC[10,9] results 
can be understood from the normal mode sketch in Fig. 3. At the MCSCF level the 
coupling of the two hydrogen rocking motions causes the two OH bond dipoles to move 
out of the plane in opposite directions, reducing the overall change in dipole moment 
and causing the intensity to drop from 128 to 42 km/mol. In the case of mode 5, the 
coupling of the rocking motions causes both bond dipoles to move out of plane in the 
same direction, increasing the overall change in dipole moment and causing a small 



intensity increase from ROHF to MC[10,9]. 
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IV. CONCLUSIONS 

Theoretical values for the infrared absorption spectrum of both the covalent and 
hydrogen bonded isomers of the complex formed between the hydroxyl and hydroperoxyl 
radicals have been presented. For the covalent HOOOH species, the vibrational 
frequencies of the stiff modes (bends and stretches) have been predicted with an 
expected accuracy of a few percent. Consideration of the spectra of possible interferents 
(Figs. 1 and 2), suggests that the relatively strong absorptions due to the bending motions 
and antisymmetric stretch of the 0-0-0 backbone may be the most easily identified and 
characterized. These are predicted to occur at 518 and 756 cm' 1 . The stronger 
absorptions due to the torsional motions of HOOOH may also be useful for this purpose, 
but it is expected that these vibrations are not as well treated within the harmonic 
approximation and may be less accurately predicted. 

Formation of the hydrogen bonded H0-H0 2 complex gives rise to five new 
interradical modes of vibration. Some of these are predicted to be quite intense and 
may be useful for detection of the species. Noteworthy in this regard are the out-of- 
plane motions of the hydrogen atoms, denoted v 2 and v 5 , which at the best 
computational level are predicted to lie in the region 700-800 cm' 1 . However, these 
frequencies show a very strong dependence on the theoretical level and are probably not 
well treated in the harmonic approximation, indicating that their predicted values should 
not be considered highly accurate. These uncertainties considerably lessen the utility of 
these modes for identification purposes. The frequency shifts of the internal vibrational 
modes, especially the HO stretch and HOO bend in the H0 2 radical, are probably more 
accurately predicted and should be most useful for detection of the HO—HOO species. 
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The large intensity increase predicted for the hydroxyl radical stretching vibration upon 
complexation may also make that mode useful. 

In addition to predicting features of the absorption spectra of these two species, 
we note two other interesting features of this study. In the case of the covalent HOOOH 
species, the magnitude and relative order of the 0-0-0 backbone symmetric and 
antisymmetric stretching mode frequencies are very sensitive to the level of theoretical 
treatment. This is quite similar to the situation previously found for the ozone molecule. 
In that case, the interpretation was complicated by the difficulty in separating effects due 
to the biradical character of ozone’s electronic ground state from other aspects of the 
electron correlation problem. In the present case, however, the closed shell singlet 
reference state is a good zeroth order wave function and similar difficulties still arise 
here. Because enlarging the active space to include the oxygen atom lone pairs has so 
much larger effect on v 5 (antisymmetric stretch) than on v 4 (symmetric stretch), it seems 
clear that this additional correlation becomes increasingly important as the molecule is 
distorted away from C 2 symmetry. An attempt was made to identify specific 
configurations that increased in importance in the CAS wave function as the molecular 
framework was distorted small amounts along the antisymmetric -OOO- stretch. 

However, examination of these Q and C 2 wave functions did not reveal any important 
changes among the configurations with | c s | ^ 0.025. Recalling that inclusion of both 
lone pairs on each oxygen atom, rather than just one, provided further improvement in 
the ozone case 46 , we anticipate that v s in the present study would also be improved by 
further expansion of the active space. However, we are not able to carry out that larger 
study at this time. 
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In the case of the hydrogen bonded species, certain of the intermolecular 
vibrations have predicted frequencies that are very sensitive to the inclusion of electron 
correlation via CASSCF, with v 2 and v 5 increasing relative to the ROHF results by 
factors of 2-3. These out-of-plane modes, which mostly involve hydrogen atom motions, 
move from quite low frequencies (100-400 cm' 1 ) at the RHF level into the 700-800 cm' 1 


region at the MC[10,0] level. 
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METHOD 

Rqh 

Rqo 

®ooo 

® HOO 

* HOOO 


DZP Basis 


RHF 

0.950 

1.369 

107.5 

103.8 

81.4 

MC[8,8] 

0.975 

1.460 

105.1 

100.0 

81.2 

MC[14,11] 

0.976 

1.457 

106.7 

100.1 

80.7 

DZP( + ) Basis 
RHF 

0.950 

1.368 

107.6 

103.9 

81.8 

MC[8,8] 

0.975 

1.458 

105.2 

100.2 

81.5 

MC[14,11] 

0.976 

1.456 

106.8 

100.4 

81.4 

DZP(2d) Basis 
RHF 

0.948 

1.368 

107.6 

103.7 

80.5 

MC[8,8]. 

0.973 

1.458 

105.1 

100.0 

79.0 

MC[14,11] 

0.974 

1.456 

106.7 

100.2 

79.8 

6-31G’/MP2 b 

0.972 

1.439 

106.3 

100.2 

78.1 


a. Units are Angstroms and degrees. The molecule has Q symmetry. R AB is the AB 
bond distance; 0 ABC is the included angle formed by the AB and BC bonds; $ HOOO 
is the dihedral angle formed by the OOO and HOO planes. 

b. Reference 22 ; see also Reference 19 for similar results obtained with the slightly 
larger 6-31G" basis set. 






TABLE II. Vibrational Frequencies (in cm' 1 ) for Covalent H 2 0 3 a . 
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METHOD 

v ~i 

V ~2 

V ~3 

Normal Mode 
v ~4 *5 v 6 

v 7 

00 

v 9 


DZP Basis 


RHF 

397 

426 

627 

1142 

1182 

1565 

1572 

4128 

4132 

MC[8,8] 

335 

363 

509 

808 

816 

1404 

1422 

3739 

3741 

MC[ 14,11] 

368 

425 

509 

833 

725 

1391 

1403 

3720 

3722 

DZP( + ) Basis 
RHF 

385 

414 

629 

1143 

1181 

1562 

1572 

4123 

4126 

MC[8,8] 

323 

348 

508 

808 

818 

1403 

1418 

3736 

3738 

MC[14,11] 

337 

391 

509 

853 

724 

1424 

1433 

3715 

3711 

DZP(2d) Basis 
RHF 

394 

433 

629 

1125 

1159 

1578 

1586 

4125 

4128 

MC[8,8] 

336 

366 

51f 

801 

814 

1419 

1429 

3737 

3739 

6-31G** Basis 
MP2 b 

366 

417 

537 

904 

829 

1395 

1399 

3792 

3796 

MP4(SDTQ) 

364 

417 

520 

868 

783 

1383 

1387 

3781 

3785 

Exptl c 



500 

855 

755 






a. Theoretical frequencies are harmonic. Experimental frequencies are observed 

fundamentals in condensate. See text for qualitative description of normal modes. 


b. Reference 19 


c. 


Reference 38 
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TABLE III. Calculated Harmonic Vibrational Frequencies for Hydrogen Peroxide (in cm" 1 ). 
Ratios to Observed Fundamentals 3 are Shown in Parentheses. DZP Basis Set used Throughout. 




Normal Mode Description 



METHOD 

OO str 

® HOO,symm 

® HOO.asym 

HO str symm 

HO str asym 

RHF 

MC[6,6] 

MC[10,8] 

MC[14,10] 

1168 (1.33) 
826 (0.94) 
833 (0.95) 
824 (0.94) 

1604 (1.16) 
1445 (1.04) 
1438 (1.04) 
1411 (1.02) 

1474 (1.16) 
1317 (1.04) 
1309 (1.03) 
1284 (1.01) 

4169 (1.16) 
3763 (1.05) 
3741 (1.04) 
3737 (1.04) 

4169 (1.15) 
3766 (1.04) 
3746 (1.04) 
3741 (1.04) 

MP2 b 

932 (1.07) 

1435 (1.03) 

1308 (1.03) 

3864 (1.07) 

3863 (1.07) 


a. Experimental Fundamentals used in determining ratios: 875, 1387, 1266, 3599, and 
3611 cm' 1 , respectively. See Reference 40. 


b. Reference 39. 
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TABLE IV. Absorption Intensities (in km/mol) a for Covalent H 2 0 3 b . 


METHOD 

v ~i 

*2 

V ~3 

DZP Basis 

RHF 

172.0 

132.2 

36.6 

MC[8,8] 

137.8 

114.5 

29.6 

MC[14,11] 

123.0 

114.9 

44.8 

DZP(+) Basis 

RHF 

176.2 

133.5 

34.2 

MC[8,8] 

142.0 

115.4 

24.9 

MC[14,11] 

126.3 

110.3 

31.3 

DZP(2d) Basis 

RHF 

153.8 

122.1 

29.6 

MC[8,8] 

120.4 

109.4 

24.1 

6-31G" 

MP2 

125.7 

134.4 

40.2 


Normal Mode 


V 4 

v 5 

v 6 

v ~7 

V 8 

v 9 

3.2 

37.4 

34.6 

56.2 

121.3 

9.0 

7.2 

15.6 

37.6 

63.4 

36.8 

2.5 

7.2 

58.3 

35.9 

57.9 

34.6 

2.5 

2.7 

39.3 

34.6 

51.1 

134.8 

8.9 

5.9 

17.3 

38.0 

58.3 

43.1 

3.0 

6.8 

61.7 

43.1 

63.0 

72.7 

10.1 

2.9 

34.2 

36.3 

60.0 

125.5 

9.7 

6.3 

14.8 

39.7 

65.5 

43.9 

3.6 

9.0 

59.7 

36.5 

63.0 

53.4 

4.3 


a. 


b. 


1 km/mol = 42.255 D 2 /amu-A 2 . 

See text for qualitative description of normal modes. 
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TABLE VI. Harmonic Vibrational Frequencies (in cm' 1 ) for 3 A’ Hydrogen-bonded HO— HOO 
Complex. Results for Uncomplexed Species Shown in Italics. 



METHOD 

*1 

*2 

*3 

Normal Mode 3 
v ~4 ^5 v 6 

v 7 

v 8 

*9 


ROHF 

DZP 

61 

148 

178 

333 

403 

1274 

1270 

1652 

1594 

4051 

4104 

4072 

4082 

DZP( + ) 

56 

121 

169 

314 

383 

1276 

1272 

1650 

1596 

4053 

4102 

4077 

4085 

DZP( + ,2d) 

60 

134 

165 

339 

387 

1247 

1242 

1664 

1611 

4036 

4095 

4063 

4079 

MCriO,91 b 

DZP( + ,2d) 

104 

723 

176 

363 

772 

1065 

1048 

1532 

1474 

3670 

3692 

3686 

3680 

MP2 

DZP(+,2d) 

138 

268 

215 

487 

467 

1262 

1208 

1533 

1454 

3527 

3708 

3752 

3828 


a. Normal mode identifications: v 134 (intermolecular in-plane (a’) modes); v 24 (intermolecular 
out-of-plane (a") modes); v 6 (0-0 stretch); v 7 (HOO bend); v 8 (HO stretch in H0 2 ); and 
v 9 (HO stretch in HO). These labels are used for similar modes in the isolated species. 


b. MC[10,9] for complex, MC[7,6] for H0 2 , and MC[3,3] for HO. 
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TABLE VII. Calculated Frequency Shifts (in cm' 1 ) for 3 A’ HO—HOO upon Complexation. 3 



METHOD 

H0 2 Modes(v 6 -v 8 ) 

OO Str HOO Bend HO Str 

OH Mode(v 9 ) 
OH Str 


ROHF 
DZP Basis 

4 

58 

-53 

-10 

DZP(+) 

4 

54 

-49 

-7 

DZP(+,2d) 

5 

53 

-59 

-16 

MCSCF 17 

DZP(+,2d) 

17 

58 

-22 

6 

MP2 

DZP(+,2d) 

54 

79 

-181 

-76 


a. Calculated as: v j(complex) - v (uncomplexed radical) 

b. MC[10,9] for complex, MC[7,6] for H0 2 , and MC[3,3] for HO. 
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TABLE VIII. Absorption Intensities for Hydrogen-bonded HO—HOO (in km/mol) a Results for 
Uncomplexed Species Shown in Italics. 


METHOD 

Vi 

v 2 

v 3 

Normal Mode 3 
v 4 v 5 

V ~6 

*7 

*8 

*9 

ROHF 










DZP 

44.4 

136.5 

25.8 

147.0 

221.8 

37.2 

63.0 

215.1 

60.4 







40.1 

56.2 

58.7 

23.7 

DZP( + ) 

31.7 

147.0 

20.7 

147.0 

210.9 

38.9 

54.1 

251.4 

58.7 







41.0 

51.6 

67.6 

30.4 

DZP(+,2d) 

36.3 

128.5 

24.5 

143.2 

186.3 

36.8 

56.6 

241.7 

66.8 







39.3 

55.4 

70.6 

34.2 

,MC[10.9] b 










DZP(+,2d) 

24.9 

41.8 

48.2 

125.1 

211.3 

23.1 

72.7 

141.6 

46.1 







28.3 

71.0 

39.3 

3.5 

MP2 










DZP( + ,2d) 

68.9 

130.2 

52.3 

162.5 

160.9 

88.7 

50.7 

305.1 

67.5 







98.7 

33.1 

55.6 

22.4 


a. Normal mode identifications: v 134 (intermolecular in-plane (a’) modes); v 2 ^(intermolecular 
out-of-plane (a") modes); v 6 (0-0 stretch); v 7 (HOO bend); v g (HO stretch in HO z ); and 

v 9 (HO stretch in HO). These labels are used for similar modes in the isolated species. 

b. MC[10,9] for complex, MC[7,6] for H0 2 , and MC[3,3] for HO. 
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Figure 1. Predicted Spectrum (200 - 1800 cm' 1 ) of Covalent HOOOH with Possible Interfering 
Species. Experimental Frequencies: H0 2 -, H 2 0 2 , and H 2 0 . 

Figure 2. Predicted Spectrum (3200 - 3800 cm' 1 ) of Covalent HOOOH with Possible Interfering 
Species. Experimental Frequencies: HO-, H0 2 % H 2 0 2 , and H 2 0 . 

Figure 3. Normal Mode No. 2 for 3 A’ HO-HOO: a) ROHF (v 2 = 134 cm' 1 ); b) MC[10,9] 

(v 2 = 723 cm' 1 ). 

Normal Mode No. 5 for 3 A’ HO-HOO: a) ROHF (v 5 = 387 cm' 1 ); b) MC[10,9] 

(v 5 = 772 cm' 1 ). 
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Figure 3 - Jackets 
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Figure 4 - Jackels 



